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A  significant  portion  of  the  power  loss  in  a  fuel  cell  stack  can  be  attributed  to  the  contact  resistance 
between  the  gas  diffusion  layer/bipolar  plate  and  the  current  collector/bipolar  plate  interfaces,  especially 
when  an  oxide  layer  is  formed  on  a  stainless  steel  bipolar  plate.  Researchers  have  already  studied  methods 
to  decrease  the  contact  resistance  between  fuel  cell  components,  but  never  has  a  theoretical  contact 
mechanics  model  been  applied  to  contact  resistance  problems  in  fuel  cells.  Therefore,  the  purpose  of  this 
research  is  to  utilize  a  theoretical  contact  mechanics  model  in  order  to  study  real  contact  area  in  fuel  cell 
components  as  a  function  of  surface  roughness,  material  properties,  and  clamping  force.  Specifically,  the 
effects  of  bipolar  plate  surface  roughness,  coating  thickness  of  the  gold  plating  on  the  current  collector, 
and  the  clamping  force  on  real  contact  area  have  been  studied.  It  was  found  that  smoother  materials, 
thicker  gold  coating,  and  higher  clamping  force  resulted  in  a  higher  real  percentage  contact  area. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

A  fuel  cell  is  an  energy  conversion  device  that  generates  electric¬ 
ity  through  the  electrochemical  reaction  of  hydrogen  and  oxygen, 
producing  only  water  as  a  byproduct.  The  decrease  of  fossil  fuel 
resources  coupled  with  the  increased  risk  of  global  warming  is 
causing  fuel  cell  technology  to  be  an  area  of  rising  interest.  Of 
all  of  the  different  types  of  fuel  cells,  proton  exchange  membrane 
(PEM)  fuel  cells  are  the  ones  that  have  attracted  the  most  inter¬ 
est,  due  to  a  low  operating  temperature,  short  start-up  time,  high 
power  density,  and  nearly  zero  emissions.  Thus,  the  PEM  fuel  cell 
would  be  desirable  in  transportation  devices  where  pollution  and 
temperature  are  major  issues. 

PEM  fuel  cells,  see  Fig.  1,  have  great  potential  to  be  environ¬ 
mentally  friendly,  but  there  are  still  several  problems  that  need  to 
be  solved  before  they  can  be  mass-produced.  The  contact  resis¬ 
tances  between  the  gas  diffusion  layer/bipolar  plate  interface  and 
the  current  collector/bipolar  plate  interface  contribute  to  a  signifi¬ 
cant  power  loss  in  a  fuel  cell  [1  ].  The  contact  resistance  is  affected 
by  the  clamping  pressure,  the  surface  geometry  and  material  prop¬ 
erties  of  contact  components. 

Several  researchers  have  studied  the  contact  resistance  of  fuel 
cell  components.  The  effect  of  gas  diffusion  layer  clamping  pres¬ 
sure  on  fuel  cell  performance  was  studied  by  Lee  et  al.  [2].  They 
suggested  that  fuel  cell  performance  can  be  affected  by  changes 


*  Corresponding  author.  Tel.:  +1  248  370  2224;  fax:  +1  248  370  4416. 
E-mail  address:  wang@oakland.edu  (X.  Wang). 


in  electrical  contact  resistance  due  to  varying  pressures.  Too  much 
pressure  can  result  in  damage  to  the  gas  diffusion  layer,  and  too 
little  pressure  can  cause  high  contact  resistance.  Determining  the 
optimal  assembly  pressure  is  crucial  in  order  to  reduce  power  loss  in 
a  fuel  cell.  Jung  et  al.  [3  ]  studied  contact  resistance  for  different  flow 
field  combinations  in  fuel  cells.  Wang  et  al.  [4]  considered  stainless 
steel  as  a  possible  material  for  a  bipolar  plate  by  testing  the  con¬ 
tact  resistance  between  stainless  steel  and  carbon  paper.  Mishra 
et  al.  [5]  compared  results  from  predicted  contact  resistance  and 
the  experimental  contact  resistance  of  different  gas  diffusion  layer 
materials  with  variable  contact  pressure.  They  utilized  a  fractal 
asperity  model  to  determine  contact  resistance  based  on  clamping 
pressure,  material  properties,  and  surface  geometry.  Ihonen  et  al. 
[6]  measured  the  contact  resistance  of  un-plated  and  plated  stain¬ 
less  steel.  Zhang  et  al.  [1]  developed  two  semi-empirical  methods 
for  estimating  the  contact  resistance  between  a  bipolar  plate  and 
gas  diffusion  layer.  They  reported  good  agreement  between  their 
predictions  and  their  experimental  contact  resistance  results. 

It  is  well  known  that  contact  resistance  depends  strongly  upon 
real  contact  area  and  there  have  been  many  successful  theoretical 
models  developed  to  determine  real  contact  area  as  a  function  of 
surface  roughness  and  material  properties  [7].  To  date,  these  mod¬ 
els  have  not  been  applied  to  contact  resistance  problems  in  fuel 
cells.  Therefore,  the  objective  of  this  research  is  to  utilize  a  pre¬ 
viously  developed  contact  mechanics  model  to  study  real  contact 
area  at  the  current  collector/bipolar  plate  interface  in  a  fuel  cell  as  a 
function  of  surface  roughness  and  material  properties  and  clamp¬ 
ing  pressure.  Upon  successful  application  of  the  contact  mechanics 
model  to  the  current  collector/bipolar  plate  interface,  the  model 
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Nomenclature 

E 

equivalent  modulus  of  elasticity 

g(x,y) 

gap  between  contact  bodies 

go 

initial  gap  between  the  reference  planes  of  two  sur¬ 
faces 

G 

equivalent  shear  modulus 

p(?.  v) 

contact  pressure  distribution 

Pj  (x,y) 

contact  pressure  for  a  given  point 

<?(?,  ri) 

tangential  force 

Ra 

average  roughness 

Rt 

roughness  total 

R, 

root  mean  square  roughness 

Rz 

Ten  point  height 

u(x,y) 

vertical  surface  displacement 

Gree/c  symbols 

8 

surface  approach 

V 

Poisson’s  ratio 

will  be  applied  to  the  gas  diffusion  layer/bipolar  plate  interface  in 
future  research  after  the  surface  roughness  distribution  for  the  gas 
diffusion  layer  becomes  available. 

2.  Contact  mechanics  model 

The  contact  mechanics  model  used  in  this  research  was  devel¬ 
oped  by  Hu,  Barber,  and  Zhu  [8].  When  two  rough  surfaces  are  in 
contact,  the  gap  between  the  deformed  surfaces  in  contact  g(x,  y) 
can  be  expressed  as 

g(x,y)  =  g0-8  +  z{x,y)  +  u{x,y)  (1) 

where  g0  is  the  initial  separation  between  two  surfaces  in  contact 
when  there  is  no  external  load.  When  a  certain  load  is  applied,  the 
two  surfaces  are  deformed,  and  they  have  an  approach  distance  of 
8.  The  surface  roughness  z(x,  y)  is  calculated  as 

z(x,y)  =  z!(x,y)  +  z2(x,y)  (2) 


Current  Collector 


Anode  Cathode 

H2-»2H++2e‘  1/202+2Ht+2e^H2° 


Current  Collector 


where  the  subscripts  1  and  2  represent  the  surfaces  1  and  2, 
respectively.  u(x,  y)  is  the  vertical  surface  displacement  due  to  a 
distributed  normal  pressure  p(§,  rj)  and  tangential  force  q(f,  rj),  and 
it  can  be  obtained  by  the  following  equation: 


u(x,y) 


P(£  P) 

\/(x-?)2  +  (y->?)2 


d|  dr, 


S-x 


\/(x-£)2+(y-J7)2 


df  dr] 


(3) 


where  E  is  the  equivalent  Young’s  modulus  and  G  is  the  equivalent 
shear  modulus  given  by 


1  l/1-uj  1— «22\ 
£  2  ^  Ei  E2  ) 


(4) 


Fig.  1.  Schematic  of  PEM  fuel  cell. 


Fig.  2.  Flowchart  for  numerical  procedure. 
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1  _  1  /(l+u1)(l-2u1)  (1  +  u2)(1-2u2)\  (5) 

G  2  V  E\  E2  ) 

where  v  is  Poisson’s  ratio. 

It  is  assumed  that  the  friction  force  is  proportional  to  the  nor¬ 
mal  force,  and  the  coefficient  of  friction,  /x,  is  a  constant.  Then,  the 
force-displacement  relation  can  be  given  by 

U{X’y)=^E  f  l  P^1l] 

_ 1 _ E  jfjl) 

V(£-x)2+(r?-y)2  2G  v'(|-x)2+(j/-y)2 

(6) 

Eq.  (6)  is  not  a  complete  solution  since  it  requires  the  knowledge  of 
either  pressure  distribution  or  displacement  distribution,  neither  of 
which  is  known  in  advance  except  the  load  and  surface  geometry. 
Therefore,  the  additional  constraints  needed  to  solve  Eq.  (6)  are 

g(x,y)>  0  p(x,y)>  0  (7) 

where  the  assumption  of  no  tension  stress  and  no  surface  penetra¬ 
tion  leads  to  a  non-negative  constraints  for  both  surface  gas  and 
contact  pressure. 

Eq.  (6)  with  constraints  (7)  is  solved  through  a  fast  Fourier  trans¬ 
formation  (FFT)  based  approach  [8,9]  and  a  conjugate  gradient 
method  (CGM)  [9,10].  The  flowchart  of  the  numerical  procedure  is 
shown  in  Fig.  2.  The  real  contact  pressure  distribution  between  the 
current  collector  and  the  bipolar  plate  can  be  obtained  through  the 
model.  Based  on  the  pressure  distribution,  the  real  area  of  contact 
is  calculated,  which  is  the  sum  of  the  areas  with  positive  pressure 
acting  on  them.  The  real  contact  ratio  is  the  ratio  between  the  real 
area  of  contact  and  the  nominal  area  of  contact. 

The  inputs  required  for  the  model  are:  Young’s  modulus  and 
Poisson’s  ratio  for  the  current  collector  and  bipolar  plate,  clamping 
load,  nominal  contact  area,  surface  roughness  data  and  coating  film 
thickness.  The  inputted  surface  roughness  data  consists  of  digital 
surface  height  data  measured  with  a  WYKO  optical  surface  profiler. 

3.  Results 

The  contact  mechanics  model  is  utilized  to  study  the  effects  of 
roughness  of  the  graphite  plate,  thickness  of  the  gold  coating  of  the 
current  collector,  and  normal  load  between  the  current  collector 


Table  1 

Material  properties  for  current  collector  and  bipolar  plate 


and  the  bipolar  plate  on  real  contact  area.  A  PEM  fuel  cell  typically 
makes  use  of  a  copper  current  collector  coated  with  gold  and  a 
graphite  bipolar  plate.  Properties  for  these  materials,  which  were 
inputted  to  the  contact  mechanics  model,  are  given  in  Table  1.  A 
clamping  load  of  44.3  N  was  applied  on  a  nominal  contact  area  of 
4.94  mm2.  This  produces  a  clamping  pressure  higher  than  that  typ¬ 
ically  used  in  fuel  cell  stacks.  However,  the  trends  seen  in  this  study 
should  be  applicable  in  any  pressure  range. 

3.1.  Effect  of  surface  roughness 

Three  different  bipolar  plate  surface  roughness  values  were 
used.  Surface  1  was  an  actual  three-dimensional  surface  roughness 
measurement  of  a  ground  graphite  bipolar  plate  from  an  actual 
PEM  fuel  cell  shown  in  Fig.  3.  To  obtain  a  range  of  roughness  values, 
surface  2  and  surface  3  height  values  were  obtained  by  multiply¬ 
ing  surface  1  values  by  a  factor  of  5  and  10,  respectively.  Surface 
1  was  found  to  have  an  arithmetic  average  roughness,  Ra,  value 
of  1.1  [xm.  Thus,  surface  2  had  a  Ra  value  of  5.5  [xm,  and  surface  3 
had  a  Ra  value  of  11.0  |xm.  The  roughness  for  the  current  collector 
was  measured  from  an  actual  current  collector  and  was  found  to 
have  a  Ra  value  of  0.073  |xm  as  shown  in  Fig.  4.  The  height  data 
from  these  roughness  measurements  were  inputted  into  the  com¬ 
puter  simulation  and  the  output  was  the  real  contact  area  between 
the  bipolar  plate  and  the  current  collector.  The  calculated  real  per¬ 
centage  contact  areas  calculated  from  the  contact  mechanics  model 
versus  surface  roughness  are  shown  in  Fig.  5.  The  calculated  areas 
are  the  sum  of  all  contact  points  between  the  bipolar  plate  and  the 
current  collector.  As  expected,  the  real  percentage  contact  area  was 
the  highest,  78.3%,  for  the  smoothest  plate  followed  by  the  inter¬ 
mediate  roughness  specimen  with  43%  contact  area,  and  finally  the 
roughest  specimen  with  a  contact  area  of  27.4%. 

3.2.  Effect  of  coating  thickness 

Three  different  coating  thicknesses  on  the  current  collector  were 
studied.  Surface  1  had  a  coating  thickness  of  7.5  [xm,  surface  2  had  a 


1.8 

1.6 
1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 

0-0  ,  i  i  .  i 

0.0  0.5  1.0  1.5  2.0  2.4 

Fig.  3.  Three-dimensional  surface  trace  of  graphite  plate.  (For  interpretation  of  the  references  to  color  in  the  artwork,  the  reader  is  referred  to  the  web  version  of  the  article.) 


Surface  Statistics: 
Ra:  1.11  um 
Rq:  1 .46  um 
Rz:  15.39  um 
Rt:  22.72  um 
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Surface  Statistics: 

Ra:  72.59  nm 
Rq:  178.54  nm 
Rz:  10.36  um 
Rt:  15.89  um 


Fig.  4.  Three-dimensional  surface  trace  of  current  collector.  (For  interpretation  of  the  references  to  color  in  the  artwork,  the  reader  is  referred  to  the  web  version  of  the 
article.) 


coating  thickness  of  12.5  |xm,  and  surface  3  had  a  coating  thickness 
of  25  |jim.  The  coating  thickness  was  inputted  into  the  computer 
model  and  the  output  was  the  real  percentage  contact  area  between 
the  current  collector  and  the  bipolar  plate  as  shown  in  Fig.  6.  The 
highest  real  contact  area,  80.2%,  was  obtained  with  a  coating  thick¬ 
ness  of  25  [xm,  followed  by  the  intermediate  coating  thickness  of 
12.5  |xm  with  a  78.8%  contact  area,  and  the  lowest  contact  area, 
72.5%,  was  obtained  with  the  7.5  |xm  coating  thickness. 

3.3.  Effect  of  clamping  force 

Two  different  clamping  forces  were  used.  One  clamping  force 
was  44.3  N  and  the  other  was  92.0  N.  These  were  inputted  into 
the  computer  simulation  and  the  output  was  the  real  contact  area 
between  the  bipolar  plate  and  the  current  collector.  As  expected, 
the  highest  real  contact  area  (83.0%)  was  obtained  with  the  higher 


Roughness  (pm) 

Fig.  5.  Real  contact  area  percentage  vs.  roughness  of  graphite  plate. 


Fig.  6.  Real  contact  area  percentage  vs.  gold  coating  thickness. 


clamping  load.  The  clamping  load  of 44.38  N  had  a  real  contact  area 
of  78.3%. 

4.  Conclusions 

In  this  paper,  a  contact  mechanics  model  was  used  to  determine 
the  effects  of  surface  roughness,  coating  thickness,  and  clamp¬ 
ing  force  on  the  real  contact  area  between  the  current  collector 
and  bipolar  plate.  From  the  data,  the  following  conclusions  are 
drawn: 

1.  When  the  surface  is  rougher,  there  is  a  lower  real  percentage 
contact  area.  It  is  important  that  the  bipolar  plate  and  current 
collector  are  finished  to  be  as  smooth  as  economically  possible. 
Also,  graphite  is  known  to  be  highly  porous,  so  graphite  bipo¬ 
lar  plate  materials  should  be  processed  and  selected  with  the 
minimum  porosity. 
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2.  The  real  contact  area  increases  with  thicker  gold  coatings  on  the 
current  collector.  This  is  likely  due  to  the  lower  elastic  modulus 
of  the  gold  coating  as  compared  to  the  copper  base  material. 
The  gold  elastically  deforms  more  than  copper  with  a  resulting 
increase  in  real  contact  area. 

3.  A  larger  contact  load  leads  to  higher  real  contact  area.  The  largest 
clamping  load  possible  should  be  used  in  PEM  fuel  cells  without 
damaging  the  gas  diffusion  layer. 

4.  The  contact  mechanics  model  was  successfully  used  to  deter¬ 
mine  the  real  contact  area  between  the  current  collector  and 
the  bipolar  plate.  This  model  has  the  potential  to  be  utilized  at 
other  contact  interfaces  in  PEM  fuel  cells  and  to  study  the  effects 
of  new  materials  and  various  contact  geometries  in  PEM  fuel 
cells. 
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